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Abstract — In this study, a novel computational framework 
for solving a class of optimal control problem for systems 
governed by non-linear Hammerstein integral equations is 
presented. A theorem for the convergence and the validity of the 
approach is also given in detail. Numerical experiments and 
comparisons with exact solutions confirm the efficiency and the 
accuracy of the proposed technique. 

Index Terms — Optimal Control problem, Non-linear 
Hammerstein integral equations, Iterative methods, 
Approximate-analytical solution. 


I. Introduction 

Iterative schemes and optimizations related to integral 
equations are two prominent fields of research in applied 
science and engineering. The major purpose of optimization 
is to determine procedures of how optimally change or 
influence real systems to achieve a desired result. This 
requires to realize large-scale optimization strategies with 
increasing complexity which in turn motivates the 
development of numerical techniques for optimization 
purposes. 

On the other hand, in mathematical formulation of physical 
phenomena, integral equations are always encountered and 
have attracted much attention. In fact, integral equations are 
as important as differential equations and appear in a variety 
of applications in many fields including continuum 
mechanics, potential theory, geophysics, electricity and 
magnetism, kinetic theory of gases, hereditary phenomena in 
biology, quantum mechanics, mathematical economics, 
population genetics, medicine, fluid mechanics, steady state 
heat conduction, and radiative heat transfer problems [1,3,6, 
7, 9, 10, 12, 13, 14, 21, 23, 25], In this way there are many 
direct and indirect numerical solution for integral equations 
such as [2, 5, 17, 17, 22], Furthermore, optimal control of 
systems governed by integral equations are momentous in 
applications such as the optimal control problem related to the 
Ornstein-Uhlenbeck process which arises from statistical 
communication theory [8]. In this paper, we focus on the 
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formulation of a class of optimal control problems governed 
by non-linear Hammerstein integral equations as follows 

Minimize J(x,u) = (1) 

subject to 

x(t) = y(t) + A, Jfc, (t, s, u(s))<t>(s, x(s))ds (2) 

+ j" o & 2 it, s, uis))Wis, x(s))ds, 

where the known function 

fo(t, xit), u(t )) e C([0,1] X R x R) and 

y(t), ®(f, xit)), Wit, x(t)),k x {t, s, uis)), k 2 {t, s,uis )) 
are given functions which are defined on the interval 
0 < t, s < 1, and can be expanded to the Taylor series about 
f = c,0 < C < 1. Besides, x(f),w(f) e C°°([0,1]) are the 
trajectory and control functions, respectively. Here, it is 
assumed that the problem (l)-(2) has a unique solution. 

Nevertheless, there is no research to solve the optimal control 
problem (l)-(2). Thus, the main purpose of this study is to 
construct an iterative scheme to obtain the approximate 
solution and also the analytical solution of the problem in the 
form of polynomial series solution. 

II. Taylor solution 

In this section, following the work of Mahmoodi [11], it is 
assumed that for a given uis) , the non-linear integral 
equation(2), for which yit), ®(f, x(f )), *+*(7, x(t)) , 
k, it,S,uis)) and k 2 it,s,uis)) are functions that have 
suitable derivatives in the interval 0 < S, t < 1 , has the 
solution x(t) in the form 
N 1 

xit) = V— x (n) ic)it-c) n , 0<t,c<l. (3) 

t^n\ 

Noting that, Eq.(3) is a Taylor polynomial of degree N at 
t = c, where X (n> ic ), n = 0,1,..., N are coefficients to be 
determined. 

Using Eq.(3), we also consider ®(7, xit)) and 
Wit, xit)) in Eq.(2) expressed in terms of Taylor 
polynomials as 
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0(t,x(t)) = Y-</>„(t-c) n , 0<f,C<l. (4) 

x V(t,x(t)) = Y- V f n (t-c) n , 0<t,c<l. (5) 

t=on\ 

Which are Taylor polynomials of degree N at t = c. The 
coefficients <f> n ,l// n ,n = are non-linear 

combinations of X (0> (c),X 0> (c),...,X (N> (c) as follows: 

K = 0(t,x(t)) 

, do „ A so 

6, = — + x (t) — 

' dt dx 

, d 2 o „ N 80 80 , u , d 2 o „ r 80 

(/) 2 =— T + x(t) — — + (x{t)) ~Y + X (0 — 
8t dx dt dx dx 

, d 3 0 . ,,,.<30(90 <3 2 0 50 

6, = —r + 2 x"(t) -+ 2(x (t)y —=-+ 

^ dt 3 dx dt dx 2 8t 

_ .... ,..<3 2 0 ,..<30<3 2 0 .. 3 5 3 $ 

3a: (t)x (t)— r + x(t)— - — r + (x(t))' —— 

dx dx dt~ dx 

+ x (f) — 

dx 


A. Matrix representation of the components 

To obtain the solution of (2) in the form of 
expression (3), we first differentiate it n times with respect 
to t : 

d n , 


x (n) (t) = y M (t) + \ -[ k x (t, s, u(s))0(s, x(s))ds 

dt n Jo 

p)t n 




d n k 2 (t,s,u(s )) 
dt n 


Here, we have 


( 8 ) 


and then analyse it as matrix representation. 

According to [15], Eq.(8) can be written as 
x <n, (t) = y M (t) + \V (n \t) + \F M {t) (9) 

where 

V (n) (t) = — f kdt,s,u(s))0(s,x(s))ds 
dt n J o 


X ¥(s,x(s))ds. 


V (0 \t ) = j\(f, s, u(s))0(s, x(s))ds, 
V (1 \t ) = ^a,5,M(y))O(f,x(0)O(y,A:(y))Jy, 


V s 0 ='¥(t,x(t)) 
d¥ „ . dT 
dt dx 

d 2 ¥ „ , d¥ d¥ . ,. » d 2 ¥ 

^ = W + X(t) ^ + (x(t)) ^ 

x d¥ 

+ x"(t)— 
dx 

a 3 T _ „. .dv d¥ _. ,. ^ d 2x v d¥ 

V / 3 = ^T + 2x (t)—— + 2(x(t)) —^ — 
dt dx dt dx dt 

n , ,, A d 2 W „ . d¥ d 2 ¥ 

+ 3x (t)x(t)— T + x(t) 

dx dx dt 

, aa3 d 3v F A d¥ 

+ ( x (t)y— T + x (0— 
dx dx 


dk x (t,s,u(s )) 

a 


dt 


dO(t x(t)) 80 

f fcjO, s, m(s)) s,u(s))x'(t) — 

dt dx 


f 


y (3 V?) = 30 g ky(t,s,u(s)) + 3 dk x (t,s,u{s )) dO(r,*(Q) 

dt dt dt 

+ 2 ^/jv dk x (t,s,u(s)) 8Q(t, x(t)) 
dt dx 


t / / xw 7^2 d 2 0 (r,x(f)) 
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f fc, (t,s,u(s)) 


d 2 ®(t,x(t )) 


dt 


•t d 3 k l (t,s,u(s )) 

» sT 
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x M (c) = y (n \c)+A±H n J m +A^_K nm y m 

m =0 m =0 

where 


®(s,x(s))ds. 


and then 

v (n \t)=Yjih l (tmt,xm (n ~ , ~' 


•* d n k l (t,s,u(s)) 
dt" 


®(s,x(s))ds. 


h .(t) = d ‘ k i (f ’ s ’ u W) \ 
dt 1 

Using Liebnitz's rule, we obtain 

V w (t) = 

[ (n-i-1) 






r‘ d n k,(t, s, u(s)) ... 

+ --- -Q>(s,x(s))ds. 

Jo dt n 

Note that, in Eq.(l 1), 

eT!;<•••)• 


Eqs.(3), (4) and (5), we obtain 

-i -1 


' ,<, V)=EE 

!=0 1=0 

{ Y 1 , pc d n k 1 (t,s,u(s)) 
"Jo A?" 


and 


w~ m ' _1) (c)4 


Uc ( s-c) m ds 


^>=e>j: 


ft d n k 2 (t,s,u(s)) 
~dt n 


Uc ( s-c) m ds. 

(13) 


Substituting Eqs.(12) and (13) in Eq.(9) gives 
x w (c) = 


-Er 


hr 


n (c) 


d"k l (t,S,u(s)) 
m! Jo dt n 

1 


l f=c (s-c)' n ds n>m 


m\ i{ 


df 


\ t=c (s-c)"'ds n<m 

(16) 


K, 


; m =—f 


1 rid n k 2 (t,s,u(s)) 


Uc ( s-c) m ds. 


(17) 


m! Jo dr" 

If we take n = 0,1,..., TV, relation (15) reduces to a system 
of TV + 1 non-linear equations for TV + 1 unknown 
coefficients x m (c), X (l) (c),...,X {N) (c) , as follows: 

X=Y + \H^ + X 1 K^ (18) 

where 


(11) 

W)' 


'y o) (c)‘ 


</>o 


Vo 


x (1) (c) 


/V) 


</>i 


V\ 

X = 

x (2) (c) 

,Y = 

y 2 ) (c) 

«D = 

A 

,t = 

Vi 


* W (c) 


y m (c) 


h 


V N 


dt" 


~Uc ( s-c)' n ds 


+ 4 E-v.,f OT 2 (M '“ W) uo-c)-rf 

dt* 

or, briefly, 


And H and K are (TV + 1 ) x (TV + 1 ) matrices 
defined in (16) and (17). The non-linear system of equations 
obtained in (18) can be solved using standard mathematics 
toolboxes as MATLAB. 


III. The iterative method and its convergence 

Let Q be the subset of the product space 
C“([0,l])xC°°([0,l]) contains all pairs (x(.),u(.)) , 
which satisfy Eq.(2). Also, let Q m n be the subset of Q 
consisting of all pairs (x m (.),U n (.j) , where u n (.) is a 
parameterized control function as the following polynomial 

U n (t) = &u (19) 

i=0 

and x m (.) is the extracted solution of the integral equation 
(2), which is considered as a polynomial of degree at most m 

xjf) = Je.( a o,«!,•••, a n )t j . (20) 

Here, e- ■ R n —> R, j = 0,1,..., m are continuous 
functions. Now, we consider the minimizing of J on Q m n 
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with {a k }l =0 as unknowns. This is obviously an 
optimization problem in n +1 dimensional space 
{(a 0 ,a l ,...,a n )eR n+1 : 

a 0 = M »(°) = U 0 = U nW = M ll 

k=0 

and J{x m ,u n ) may be considered as a function 

Suppose, be the solution of 

minimizing J on Q mn ,m = 1,2,_ ; n = 1,2,..., then the 

polynomial form of U n (.),n = 1,2,... (19) and using (20) 
allow us to apply the proposed method (in Section 2) for 
extracting polynomial solution of (20), which results in 
obtaining a sequence of trajectory functions } m=1 as 

Taylor series, and finally to achieve a minimizing sequence 
{(<(•),<(•) ■ 


approximate solutions will be compared with the exact 
solutions. 

Example 1. Let us first consider the optimal control problem 
governed by Fredholm-Hammerstein integral equation as 
follows 


Minimize J = £(x(f) -1 -t 2 ) 2 + ( u(t ) -1 - t) 2 dt, 


subject to 

x(t) = y(t ) + jV x u(s ) x 'PCs, x(s))ds. 


( 21 ) 

( 22 ) 


where y(t) = - + t + 1 and W(t, x(t)) = t + x(t) . 

The exact optimal solutions of (21)-(22) are 
x*(t) = \ + t 2 and u*(t) = l + t, 
with the optimal criterion 

J = J(x*(t), u* (t)) = 0. 

In the first iteration, i.e. for n = 1 and m = 1, we have 


{a m n }“ n=l is a convergent sequence. 

Proof. Proof in [16, 19]. 

Now, it can be concluded that {cc mn } is 
convergent, because it is a nondecreasing and bounded from 
below sequence. 

Theorem 1. If lim« m , n = a then a = inf q J ■ 
Proof. Proof in [16, 19]. 


m,n = 1,2,..., then 

W)" 


T 


'0 o' 

X = 

x (V] (c) 

,Y = 

0 

,K = 

0 0 


* w (c) 

l + x (1) (c) 


In the next iteration and for n = 1 and m = 2, we 


The above discussion and results can be 
summarized in a numerical algorithm for obtaining the 
approximate solutions for the optimal control (1) subject to 
Eq.(2). 

Algorithm 1. Choose £ l > 0 and £ 2 >0 for the accuracy 
of the solution. 

• Step 1. Let m, n, k = 1, u x (t) = a 0 +a 1 t , 
x x (t) = e 0 + e x (t) and a x = J(x(\),u(\)), where 

e o = e 0 (« 0 ’ a i) and = eMo’CO ■ 

• Step 2. Let m = m +1 and k = k +1 and find 
a k = inf Q m n J ■ 

• Step 3. If I OC k — Ot k | l< £ x then go to Step 4, otherwise 
go to Step 2. 

• Step 4. Let n = n +1 and k = k + 1, find 

a k = inf q J and go to Step 5. 

• Step 5. If I a k — a k _ x \< then stop, otherwise go to 
Step 4. 


IV. Numerical experiments 

Now, we show the efficiency of the method 
described using the following examples. In all examples, the 



V 0) (c)' 


1 

X = 

* (1) (c) 

x (2) (c) 

, Y = 

0 

-3.8333 

1 

0 

0 

0 

0 

0 

0 


2 a 0 + a ] a 0 + (2a,)/3 a 0 /3 + a/4 


x (0) (c) 


x (2) (c) 


So by using the iteration Algorithm 1, the numerical 
results are illustrated in Table 1 and Figures 1-4. 


Table 1: The Approximate-Analytical results for Example 1 


Iter 

n 

m 

x(t) 

u(t) 

J 

1 

1 

1 

1 

1 +t 

0.2 

2 

1 

2 

1 +t 2 

1 +t 

0 


www.erpublication.org 


























International Journal of Engineering and Technical Research (IJETR) 
ISSN: 2321-0869, Volume-3, Issue-2, February 2015 



Figure 1: Exact and approximate trajectory functions for 
Example 1, n = 1, m = 1. 



Figure 2: Exact and approximate control functions for 
Example 1, n = 1, m = 1. 



Figure 3: Exact and approximate trajectory functions for 
Example 1, n = 1, m = 2. 



Figure 4: Exact and approximate control functions for 
Example 1, n = 1, m = 2. 


Example 2. We consider the optimal control problem 
governed by Volterra-Hammerstein integral equation as 
follows 

Minimize J = f ( x(t ) - £ - £ 2 ) 2 + (w(f) -1 - t 2 ) 2 dt , (23) 

Jo 

subject to 

x(?) = _y(f) + J t x m(s) x 0(y, x(s))ds. (24) 

, , 2I . - 2 1 2 

y(t) = t-t 


®(t,x(t)) = l + t+x(t). 

The exact optimal solutions of the problem (23)-(24) 
are 

x*(t) = t + t 2 and u*(t) = l + t 2 , 

with the optimal criterion 

J = J(x*(t),u(t)) = 0. 

Using the iteration Algorithm 1, the numerical results are 
illustrated in Table 2 and Figures 5-10. 


Table 2: The Approximate-Analytical results for Example 2 


Iter 

n 

m 

x(t) 

u(t) 

J 

1 

1 

1 

t 

0.8334 + t 

0.2056 

2 

1 

2 

t + 0.9074r 2 

0.9074 + 0.8889? 

0.0086 

3 

2 

2 

t + t 2 

1 + ? 2 

0 



Figure 5: Exact and approximate trajectory functions for 
Example 2, n = 1, m = 1 . 



Figure 6: Exact and approximate control functions for 
Example 2, n = 1, m = 1. 
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Figure 7: Exact and approximate trajectory functions for 
Example 2, n = 1, m = 2. 



Figure 8: Exact and approximate control functions for 
Example 2, n = 1, m = 2. 



Figure 9: Exact and approximate trajectory functions for 
Example 2, n — 2,m — 2. 



Figure 10: Exact and approximate control functions for 
Example2, n = 2,m = 2. 


Example 3. Finally, Let us consider the optimal control 
problem governed by Volterra-Fredholm-Hammerstein 
integral equation as follows 

Minimize J = f ( x(t ) -t-t 2 ) 2 + (u(t)- 1 - t) 2 dt ,(25) 

Jo 

subject to 

x(f) = y(f) + x m(s) x <D(y, x(s))ds (26) 

+ x u(s ) x 'P(s,x(s))£/y. 

, N 15e f 2 re'fr 3 +4r 2 + 6r+ 4) 

where y(f) = t - + t - x - 1 , 

4 4 

0(f, x(f)) = 1 +1 + x(f) and W(t, x(t)) = 1 +t + x(t). 
The exact optimal solutions of (25)-(26) are 
x*(t) = t + t 2 , and u(t) = l + t, 
with the optimal criterion 

J = J(x\t),u\t)) = 0. 

Table 3 and Figures 11-14 give the numerical results for 
example 10. 



Figure 11: Exact and approximate trajectory functions for 
Example 3, n = 1, m = 1. 



Figure 12: Exact and approximate control functions for 
Example 3, n = 1, m = 1. 
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Figure 13: Exact and approximate trajectory functions for 
Example3, n = l,m = 2. 



Figure 14: Exact and approximate control functions for 
Example 3, n = 1, m = 2. 


V. Conclusion 

Optimal control problems governed by non-linear 
Hammerstein integral equations are usually difficult to solve 
analytically and so it is necessary to obtain the approximate 
solutions. The present method is effective for cases where the 
known functions have sufficient derivatives in the given 
interval. One of the advantages of this method is that the 
optimal solutions, i.e. the trajectory and control functions, are 
expressed as a Taylor series truncated at t = C . Therefore 
x(t ) and u(t) can easily be evaluated for arbitrary values of 
t with a low computation at effort. This method will not work 
for cases where the given functions do not have enough 
derivatives. An interesting feature of this method is that we 
obtain analytical solution in many cases, as shown in the 
examples. 

References 

[1] M.A. Abdou. “On the solution of linear and nonlinear integral 
equation”. Appl. Math. Comput., 146, pp. 857-871. 

[2] S. Anderyance and M. Hadizadeh, “The Computational Efficiency Of 
Walsh Approximation For Two-Dimensional Volterra integral 
equations”. Asian-European Journal of Mathematics, 4, pp. 263-270. 

[3] E. Babolian, J. Biazar and A.R. Vahidi, “The decomposition method 
applied to systems of Fredholm integral equations of the second kind”. 
Appl. Math. Comput., 148, pp. 443-452. 

[4] C. Homescu and I.M. Navon, “Optimal control of flow with 
discountinuos”. J. of computational physics, 187, pp. 660-682. 

[5] M. I. Gil, “Positivity of Green's Functions to Volterra Integral and 
Higher Order Integro-Differential Equations”. Analysis and 
Applications, 7, pp. 405-418. 

[6] S.C. Huang and R.P. Shaw, “The Trefftz method as an integral 
equation”. Adv. Eng. Software, 24, pp. 57-63. 


[7] S. Jiang and V. Rokhlin, “Second kind integral equations for the 
classical potential theory on open surface II”. J. Comput. Phys., 195, 

pp. 1-16. 

[8] T. Kailat, “Some integral equations with "Nonrational" kernels”. 
IEEE, Trans. Inf. Th. IT., 12, pp. 442-447. 

[9] P.K. Kythe and P. Puri, Computational Methods of Linear Integral 
Equations. Birkhauser Bosten, c/o Springer-Verlag, New York (2002). 

[10] D. Liang and B. Zhang, “Numerical analysis of graded mesh methods 
for a class of second kind integral equations on real line ”. J. Math. 
Anal. Appl., 294, pp. 482-502. 

[11] Y. Mahmoudi, “Taylor polynomial solution of non-linear 
Volterra-Fredholm integral equation”. Int. J. of Computer 
Mathematics, 82, pp. 881-887. 

[12] K. Maleknejad and M. Karami, “Using the WPG method for solving 
integral equations of the second kind”. Appl. Math. Comput., 166, pp. 
123-130. 

[13] K. Maleknejad and M. Karami, “Numerical solution of non-linear 
Fredholm integral equations by using multi-wavelets in the 
Petrov-Galerkin method”. Appl. Math. Comput., 168, pp. 102-110. 

[14] K. Maleknejad and Y. Mahmoudi, “Numerical solution of linear 
Fredholm integral equations by using hybrid Taylor and Block-Pulse 
functions”. Appl. Math. Comput., 149, pp. 799-806. 

[15] K. Maleknejad and Y. Mahmoudi, “Taylor polynomial solution of 
high-order non-linearVolterra-Fredholm integro-differential 
equations”. Appl. Math. Comput., 145, pp. 641-653. 

[16] O. S. Fard and M. Sanchooli,” Taylor solver for Fredholm optimal 
control problems ”. J. of Adv. Research in Differential Equations, 2, 

pp. 10-20. 

[17] J. Rashidinia and A. Parsa, “Semi-Orthogonal Spline Scalling 
Functions For Solving Hammerstein Integral Equations”. Int. J. of 
Wavelets. Multiresolution and Information Processing (IJWMIP), 9, 
pp. 427-443. 

[18] J. Saberi-Nadjafi and M. Tamamgar, “Modified Homotopy 
Perturbation Method for Solving Integral Equations ”. Int. J. of Modern 
Physics B (IJMPB), 24, pp. 4741-4746. 

[19] M. Sanchooli and O. S. Fard, “Numerical Scheme for Fredholm 
Integral Equations Optimal Control Problems via Bernstein 
Polynomials”. Australian J. of Basic and Applied Sciences. 4, pp. 
5675-5682. 

[20] I. Y.M. Stems, KJ.E. Versyck and J.F.M. Van Impe, “Optimal control 
theory: a generic tool for identification and control of (Bio-)chemical 
reactors”. Annul rewievs of control, 26, pp. 57-73. 

[21] W. Wang, “A new mechanical algorithm for solving the second kind of 
Fredholm integral equation”. Appl. Math. Comput., 172, pp. 946-962. 

[22] S. Yalcinbas and K. Erdem, “Approximate solutions of Nonlionear 
Volterra Integral Equation Systems”. Int. J. of Modem Physics B 
(IJMPB), 24, pp. 6235-6258. 

[23] S.C. Yang, “An investigation into integral equation methods involving 
nearly singular kernels for acoustic scattering”. J. Sound Vib. 234, pp. 
225-239. 

[24] T. Zdenek, T. Vosegaard, C. Kehlet, N. Khaneja, S. J. Glaser and N. C. 
Nielsen, “Optimal control in NMR spectroscopy: Numerical 
implementation in SIMPSON”. Journal of Magnetic Resonance, 197, 
pp. 120-134. 

[25] S. Zhang, Integral Equations. Chongqing press, Chongqing, China 
1987. 

Reza Mirbaqeri Jam received the MSc degree in Mathematics at 
AmirKabri University of Tehran. 

His research interests are in the areas of applied mathematics and He is 
the faculty of Islamic Azad University, Iranshahr, Iran. 

Mahmood Sanchooli is PhD student of Optimal Control at Payam Nour 
University of Mashhad, Iran. His main research interests are: Optimal 
Control, Operating Research and Fuzzy. 

Omid Solaymani Fard is associate professor in Applied Mathematics, 
and the faculty of Damqan University, Iran. His main research interests are: 
Optimal Control, Operating Research and Fuzzy. 


89 


www.erpublication.org 








